Antiferromagnetic fluctuations and d-wave superconductivity in electron-doped 

high-temperature superconductors 
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We show that, at weak to intermediate coupling, antiferromagnetic fluctuations enhance d-wave 
pairing correlations until, as one moves closer to half-filling, the antiferromagnetically-induced pseu- 
dogap begins to suppress the tendency to superconductivity. The accuracy of our approach is gauged 
by detailed comparisons with Quantum Monte Carlo simulations. The negative pressure dependence 
of Tc and the existence of photoemission hot spots in electron-doped cuprate superconductors find 
their natural explanation within this approach. 



For almost two decades, the mechanism for high tem- 
perature superconductivity has been one of the main is- 
sues in condensed matter physics. Despite an extensive 
body of theoretical work, there is at present no con- 
sensus about the mechanism. This is mainly due to 
lack of a reliable theoretical tool for a strong coupling 
problem where the value of the on-site Coulomb inter- 
action U is of the order of, or larger than, the band- 
width fH. The situation, however, appears more promis- 
ing for electron-doped high-temperature superconductors 
(e-HTSC) in which the charge gap at half-filling is 25% 
smaller than that of hole-doped cupratcs (h-HTSC), sug- 
gesting a smaller value of t/ 2]. This offers an opportu- 
nity for theories of rf-wave superconductivity based on 
weak- to intermediate-coupling approaches |2l U, y| - 

In this communication, we show that improved the- 
oretical calculations can indeed describe several aspects 
of e-HTSC that were unexplained by previous calcula- 
tions. For example, we show that the negative pressure 
derivative of the superconducting transition temperature 
Tc of e-HTSC, which contrasts with the positive pressure 
derivative of h-HTSC, can be explained. In addition, the 
hot spots observed in Angle Resolved Photoemission Ex- 
periments (ARPES) also come out of the calculation. We 
also discuss how a decrease in T,. in the underdoped re- 
gion can occur when a large pseudogap is produced by 
antiferromagnetic (AFM) fluctuations. Previous calcula- 
tions 6] predicted that would increase monotonically 
as one approaches half-filling. 

Let us first consider the results of numerical calcula- 
tions concerning the possibility of d-wave superconduc- 
tivity in the Hubbard model. In Fig.l we present a rather 
detailed survey of the dj.2_y2-wave susceptibility Xd ob- 
tained from Quantum Monte Carlo (QMC) calculations 
llBlili^ for the Hubbard model. The Hubbard model 
is characterized, as usual, by nearest-neighbor hopping 
t and on-site repulsion U. By contrast with variational 
methods, the QMC calculations are unbiased. They also 
can be performed on much larger lattices than exact diag- 
onalizations. QMC is essentially exact, within statistical 
error bars that, in Fig. 1, are generally smaller than the 
symbol size. As usual, the d^-2_y-2-w&Ye susceptibility is 
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FIG. 1: (Color) The d^2_y2 susceptibility obtained from QMC 
simulations as a function of doping and of temperature for a 
6x6 lattice. Various values of U correspond to different 
colors. The size dependence of the results is small at these 
temperatures. The Trotter step size is At = 1/10 while the 
number of measurements at each point in parameter space is 
around 10^. Measurements are grouped in blocs of 125 and 
stabilized every five steps along the imaginary-time axis. 



defined by Xd = Jo dr {TrA (t) At) with the a- wave or- 
der parameter — 9 il) '^It'^l+^x '^^^^ 1 
being over nearest-neighbors, with g (7) ±1/2 depend- 
ing on whether 7 is a neighbor along the x or the y axis. 
From now on, we work in units where fc^ = 1, ^ = 1, 
lattice spacing and hopping t are unity. The results are 
shown for various temperatures T — dopings 5 and 
interaction strengths U (shown by the various colors). 
The data clearly shows that the dome shape dependence 
of Xd is present not only for strong couphng {U > 8), but 
also at weak to intermediate coupling ([/ = 4). For weak 
coupling the dome shape occurs at temperatures that are 
sufficiently low {(3 = 4) for AFM (or spin-density-wave) 
correlations to build up. It has been known for a long 
time that these results, obtained by a numerical method 
of choice, by themselves do not suffice to decide whether 
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there is a d-wave superconducting phase in the Hubbard 
model. Indeed, the susceptibility should diverge if there 
is a phase transition. Also, at P — A the non-interacting 
model, U = 0, has a larger susceptibility than the U 
model, a fact that does not encourage optimism. 

To conclusively verify whether d-wave superconductiv- 
ity exists in this model at weak to intermediate coupling, 
one needs to reach temperatures that are an order of 
magnitude smaller than those shown in Fig.l. As is 
well known, the so-called sign-problem renders impossi- 
ble simulations at these low temperatures. To reach such 
temperatures, we use the Two-Particle Self-Consistent 
approach (TPSC) and extend it to compute su- 

perconducting correlations. The accuracy of the method 



In the above expression, Eq.(^, gd (k) is the form factor 
for the gap symmetry, while k and k' stand for both wave- 
vector and Matsubara frequencies k = (k, (2n -I- 1) ttT) 
on a square-lattice with N sites at temperature T. The 
spin and charge susceptibilities take the form x7p {q) = 
Xo(g)"^ - ^ and x~h (l) = Xaiq)~^ + ^ with xo com- 
puted with the Green function G'i-^'' that contains the self- 
energy whose functional differentiation gave the spin and 
charge vertices. The values of Ugp, Uch and (nini) are 
obtained from Ugp = U {n-^ni) / ((n^) (n|)) and from 
the local- moment sum-rule. In the pseudogap regime, 
one cannot use Ugp — U {n-^n^) / ((n-f) (n^)). Instead 
one uses the local- moment sum rule with the zero 
temperature value of (jt-i^-j,) obtained by the method of 
Ref. 01 that agrees very well with QMC calculations 

(2) 

at all values of U. Also, Go- contains self-energy effects 
coming from spin and charge fluctuations, as described 
in previous work 0, . 

The effective interaction in the particle-particle chan- 
nel mediated by AFM fluctuations is represented by the 
second term of Eq.(^. It becomes sizeable only after spin 
fluctuations have become large. Eq.(^) thus contains two 
leading effects, namely spin and charge fluctuations in- 
fluence the magnitude of the effective interactions in the 
particle-particle channel and they also decrease the life- 

(2) 

time of particles that pair (through Ga )■ The latter 
effect is generally detrimental to superconductivity while 
the former may favor pairing. 

The explicit expression for Xd, Eq.JQ), allows us to find 



for spin fluctuations and self-energy has already been 
proven 0, by comparisons with QMC data. In par- 
ticular, there is a pseudogap of AFM origin at a crossover 
temperature Tx- 

Briefly speaking, to extend TPSC to compute pair- 
ing susceptibility, we begin from the Schwinger-Martin- 
Kadanoff-Baym formalism with both diagonal 0, 0| 
and off-diagonal 0| source flelds. The self-energy is ex- 
pressed in terms of spin and charge fluctuations and the 
irreducible vertex entering the Bethe-Salpeter equation 
for the pairing susceptibility is obtained from functional 
differentiation. The final expression for the c?-wave sus- 
ceptibility is. 
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analytically which gap symmetry is enhanced or sup- 
pressed by AFM fluctuations. Indeed, since near half- 
filling AFM fluctuations are strongly peaked at k' — 
k = Q (commensurate or incommensurate), the sign of 
/ = —gd (k -|- Q) / gd (k) and the magnitude of gd (k) near 
the Fermi wave vector ki? determine the most favorable 
gap symmetry. Within a spin-singlet subspace, s-wave 
and dxiy-wave symmetries are suppressed since / < 0. 
Extended s-wave symmetry has / > 0, just like (i^2_j,2- 
wave, but its form factor is much smaller near k^, so we 
take gd (k) — (cosfc^ — cosfcj^). 

Let us first verify the accuracy of this approach by 
comparing, in Fig. 2, the QMC results for Xd, shown by 
symbols, with those of the generalized TPSC approach, 
Eq.(P), indicated by the sohd fine. The case U = 0, 
/3 4 is for reference. Fig. 2 demonstrates that the 
approach, Eq.QJ, agrees very well with QMC results for 
Xd U — A. The agreement improves for lower values 
of U . When the interaction strength reaches the inter- 
mediate coupling regime, J7 = 6, deviations of the order 
of 20 to 30% may occur but the qualitative dependence 
on temperature and doping remains accurate. The inset 
shows that previous spin-fiuctuation calculations (FLEX) 
in two dimensions ^ ^ deviate both qualitatively and 
quantitatively from the QMC results. More specifically, 
in the FLEX approach Xd does not show a pronounced 
maximum at finite doping. Moreover, it is known from 
previous work that FLEX does not show a pseudogap in 
the single-particle spectral weight at the Fermi surface 
[l3| . In TPSC the pseudogap is the key ingredient that 
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FIG. 2: (Color online) Comparisons between the d^2_y2 sus- 
ceptibility obtained from QMC simulations (see previous fig- 
ure) and from the approach described in the present work. 
QMC error bars are smaller than the symbols. Analytical 
results are joined by solid lines. Both calculations are for 
?7 = 4, a 6 X 6 lattice and four different temperatures. The 
case U = Q, 13 = A \s shown for reference. The size depen- 
dence of the results is small at these temperatures. The inset 
compares QMC and FLEX at C/ = 4, /3 = 4. 



leads to a decrease in in the underdoped regime. 

In TPSC we can understand why, as mentioned above, 
Xd is smaller than the non-interacting value in this tem- 
perature range. Indeed, the main contribution is from 
the first term in Eq. |^ which represents a pair of prop- 
agating particles that do not interact with each other. 
The contribution of the second term, which represents 
interaction through the exchange of spin and charge fluc- 
tuations, is, for (3 — 4:, about 1% at 5 = 0.5, growing to 
only 22% at (5 = 0. Hence, in this temperature range, Xd 
is smaller than the non-interacting value because of the 
decrease in spectral weight ed, uj = brought about by 
AFM self-energy effects. 

While it is impossible to do QMC calculations at lower 
T, the analytical formula for Xd, Eq.(^, can be extended 
to low T and to 256 x 256 lattice size using renormal- 
ization group acceleration and Fast Fourier Trans- 
forms. This allows us to verify whether there is d-wave 
superconductivity (d — SC) in the Hubbard model at 
weak to intermediate coupling. The complete Bethe- 
Salpeter equation would contain the possibility of repeat- 
edly exchanging spin fluctuations. Eq.Q contains only 
the first two terms, namely the zero and the one spin- 
and charge fluctuation exchange. As for the expansion 
{1 ~ x) ^ H- a;, the divergence should occur when the 
first two terms have the same magnitude. We can thus 
estimate Tc for d — SC. As usual, the Tc obtained from 
the divergence of the infinite series (Thouless criterion) 
should give an upper bound to the Kosterlitz-Thouless 
transition temperature Tkt expected in d = 2. In Fig. 3 
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FIG. 3: (Color online) Part (a) shows the contributions from 
the first term (DOS) and second term (vertex) of Eq.Q. In 
(b), our estimate of Tc using the Thouless criterion for (7 = 4 
and [7 = 6. 



(a) the first (DOS) and second (Vertex) contributions in 
Eq.(P) are plotted for C/ = 4 at /3 = 64 as a function of 
doping. The vertex part becomes larger than the first 
part over a range of 6. This signals, according to our cri- 
terion, that 0.07 < (5 < 0.13 is below at /3 = 64. Note 
that it is because the vertex part decreases much faster 
than the DOS part near half-filling that the d — SC stops 
close to half-filling, leading to a dome shape in Tc- The 
fast decrease of the vertex part near half-filling is because 
it has its strongest contribution near the Fermi surface 
where the pseudogap effect is most pronounced. 

The solid lines with symbols in Fig. 3(b) give the value 
of Tc estimated for two values of U in the intermediate 
coupling regime. The U — Q results should be viewed 
as giving the qualitative trend with increasing U. As is 
clear by now, the decrease of Tc near half-filling is caused 
by the same AFM fluctuations that enhance d — SC at 
large doping, d — SC fluctuations in our approach are 
important only between Tc and TKT^ by contrast with 
phase fluctuation theories at strong coupling jl^. Our 
results also contrast with theories where the decrease of 
Tc is driven by hidden competing broken symmetry |l9l |. 

To make more detailed connection with experimental 
results on e-HTSC, one should add second- neighbor t' 
and third-neighbor t" hopping to the Hamiltonian, as 
suggested by band-structure calculations and by ARPES. 
We perform the usual particle-hole transformation that 
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FIG. 4: (Color) Fermi surface plots obtained from energy dis- 
persion curves integrated from —0.2 to 0.1 with \t'\ — 0.175, 
\t"\ = 0.05. On both plots, 5 = 0.15, (3 = 40, system size 
128 X 128. In (a), U = 5.75 (b) U = 6.25. 

maps electron-doping of the negative t' model to hole 
doping with positive t' . As t' and t" increase, AFM fluc- 
tuations are frustrated, so pseudogap effects become less 
important and the fall of Tc on the underdoped side be- 
comes less and less pronounced. Including AFM coupling 
in the third dimension would lead to a real AFM tran- 
sition that would eventually overcome the pairing insta- 
bility. The more significant result we want to draw at- 
tention to is that, assuming that applying pressure only 
increases t (and thus decreases U/t), the data of Fig. 3 
shows that d\nTc/dP < in the weak to intermediate 
coupling regime described by our approach. This remains 
true with finite t' and t" and agrees with the experimental 
neg ative pressure dependence of Tc in these compounds 
|20|. By contrast, in h-HTSC d\nTc/dP has the opposite 
sign. If antiferromagnetism plays a role in the supercon- 
ductivity of both e-HTSC and h-HTSC, then the positive 
sign of dlnTc/ dP in the latter may be understood from 
the fact that they are in the strong-coupling regime where 
J = At'^/U increases with pressure. 

Fig. 4 shows Fermi surface maps for all wave vectors k 
in the first quadrant of the Brillouin zone. The maps are 
obtained, as in ARPES experiments on NCCO 21], from 
the integral of the single-particle spectral weight A (k, us) 
times the Fermi function over a frequency range running 
from —0.2 to -1-0.1. For U = 5.75, shown on the left-hand 
side, two hot spots are clearly apparent at the intersec- 
tion of the Fermi surface with the AFM zone boundary, 
as observed experimentally at optimal doping. The AFM 
correlation length ^ is 12 lattice spacings for this plot and 
the spin susceptibility at (tt, tt) is much larger than the 
non-interacting value. At this P and for this value of 
U, a pseudogap is observed only at the hot spots. They 
appear because the strong low-energy AFM fluctuations 
can scatter excitations at these points to other points on 
the Fermi surface separated by (tt, tt) If U is not 

large enough, there is only a decrease of spectral weight 
at the hot spots instead of a real pseudogap. By contrast, 
the right-hand side of Fig. 4 shows that if the interaction 



is too large, U = 6.25 (^ = 18), the AFM fluctuations 
scatter so strongly that a pseudogap appears everywhere 
along an arc on the Fermi surface. This confirms our con- 
tention that U cannot be too large near optimal doping in 
e-HTSC to explain the experimental results. The value 
of U, however, does have to increase with decreasing S 
so as to recover the Fermi maps observed at 5 = 10% as 
well as the Mott insulator at half- filling [illiil. 

In summary, in e-HTSC the symmetry of the super- 
conducting order parameter, the dependence of Tc on 
pressure, as well as the hot spots observed by ARPES 
at optimal doping can all be explained by the Hubbard 
model at weak to intermediate coupling. Generally, an- 
tiferromagnetic fluctuations help superconductivity until 
they are so strong that they open up a pseudogap that 
hinders d — SC. 
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